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' A giant level shift, resulted from the interaction of an electron in a spherical quantum dot with 

I zero-point oscillations of confined modes of the electric field, is divulged. The energy correction 

^SJ ■ depends on the dot radius. This size scaling of the depolarisation effect is computed semiclassically. 

A change of the optical properties of the matrix surrounding the dot provides a method to study 
^ , the shift experimentally. 
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I. INTRODUCTION 



A complete quantum dot (QD) theory, taking into account all the sophisticated physics of the object, is still a 
O [ challenge for a theorist. The main reason is that the scale of the calculation is much larger than atomic one (that 
complicates ab initio techniques) . The same time the number of particles is small to use solid state approximations 
in full extent. For example, an one-electron picture of a quantum confinement potential, arising from the conduction 
band discontinuity on the QD boundary, does not always yield accurate electron levels. 

In the paper we put forward a model to inspect the electrodynamical correction to the one-electron energy in 
a spherical QD. It was shown |^ that the similar correction occurs to be significant for a "natural quantum dot" 
Ceo- A depolarisation level shift due to the interaction with an electromagnetic field is not negligible, as it might be 
^ . thought, when taking into account localized electromagnetic modes. We present the scaling analysis of some different 

' mechanisms for the level shift (LS) and propose a possible experimental manifestation of the depolarisation effect. 
^ In order to appraise the LS a simple spherical QD model in frame of an effective mass approximation was applied. 

O How is our result sensitive to the model used? The size scaling of the depolarisation shift preserves, being dependent 
mainly on a corresponding density of states of the field, while the prefactor might be smaller within other approach, 
though it is not easy to evaluate explicitly. 
(N '. The group of full rotations, SO (3), was chosen to label the one-electron states. It is possible to perform an analytic 
quantum-mechanical calculation of the RPA response within the spherical model |^ . The massive peak of a collective 
excitation is known to show up in the spectrum, resulting from the fast coherent oscillation of the total electron 
I density of the valence states. Thus, within our model the electron-electron interaction is dealt with selfconsistently. 
CNl Of course, the number of valence electrons involved in the collective motion has not to be small. It is believed to 
fulfill for the typical QD possessing some hundreds of atoms and even more. 

A surface charge density oscillation can be thought as a confined electric field mode or a multipole surface plasmon. 
We will reflect on the shift of the electron level in the field of zero-point oscillations of the modes connected with the 
"j^ QD, which depolarisation effect is billion times stronger than of the free field zero-point oscillations, so a name "giant 
LS" is admitted. 

I The classical description of the electromagnetic surface modes, via the dielectric function of the matrix and QD 

' O material, gives the true plasmon state frequencies and will be used below. Once more, the final result does not depend 
^ too much on the computation approach. Our model sketches out the (many-body) depolarisation semiclassically, 
avoiding a lot of the routine computational intricacy. 

The paper proceeds as follows: the brief model description is next to the introduction. Then, the model will be 
applied to a 3D-plasmon as well as a free field mode, that will explain the calculation technique. However, the LSs 

deals with the confined modes those 



from these modes are too small to have an experimental importance. Section [TV 
result in much larger depolarisation 
given in respect to a possible exper: 
frequencies is allocated in Appendix 



?H result in much larger depolarisation shift. The numerical estimations and the scaling properties of the LS will be 
. . . given in respect to a possible experiment. A brief summary will follow. The calculation of the QD confined mode 



II. SEMICLASSICAL THEORY FOR ENERGY LEVEL SHIFT 

We have considered semiclassically the LS for an arbitrary shell object in The method follows the one proposed 
by Migdal Q to calculate the Lamb shift for a hydrogen-like atom. The frequency of the zero-point oscillations of 
the external field is much higher than the inverse period of the electron orbit LOp 3> tt/t. Therefore, the adiabatic 
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approximation has to be used and one divides the fast (field) and slow (electron) variables. An electron is subjected to 
short fast deflections from its original orbit in the high-frequency field of the electromagnetic wave of the zero-point 
oscillation. Then the energy shift is given by the second order perturbation theory as 

6E = {H{r + 6)- H{r)) = (vH ■ S + ^V^H 5-6+..}^, (1) 

where H{r) is the unperturbed Hamiltonian and H{r + 5) is the Hamiltonian with account for the random electron 
deflection S. The angle brackets represent the quantum mechanical average over the fast variables of the field (or, the 
same, over the random electron deflections). The perturbed Hamiltonian is expanded in series on the S and a first 
nonzero contribution is taken. 

The simplest QD Hamiltonian is considered to have only the rotational correction which is given by: 

where R is about the spherical QD radius; m is the electron mass which is supposed to be constant within the dot; 
L is the angular momentum operator. On the averaging, the first-order term disappears. So far the LS dependence 
on the QD size includes factor besides some power hidden in the mean square deflection S'^. The strength of the 
electrodynamical interaction changes with this quantity almost exactly. We will show that the dependence of S"^ in 
R is different for different electric modes (confined and free field). The giant defiection is representative for the giant 
LS and, therefore, the function 6'^{R) will be studied specifically. 



III. BULK PLASMON CONTRIBUTION TO LS 



First we consider the bulk 3D-plasmon modes that could shift the electron level. Nearly self-evidently the bulk 
plasmon shift is negligible. The mean square deflection, caused by the 3D mode (which is not confined at all), decreases 
with the QD size too rapidly. The small factor, contained in the 3D LS, comes essentially from the expression for 6'^ 
which scales as where A'' is the number of atoms in the QD. It will be explained in this section. 

Within the semiclassical approach, the deflection of the electron can be computed with the use of the Newton law: 

md^S = e£, (3) 

here e is the electron charge, £ is the field strength due to the zero-point oscillation of some mode, m is the electron 
effective mass. The square of the deflection — e^/(2m)^ J d^k Sl/oj"^ is proportional to the mean square of the 
electric field strength. The dimension of the field, D, equals 3. The field strength, in turn, can be rewritten as the 
zero-point oscillation frequency f| = 2'Khujk through the quantized field normalisation. 

The scale of the energy is given by the 3D plasmon frequency ujp = \J A.-Ke^n^x^lm. Note that the 3D plasmon 
frequency does not depend on the quantum number k and passes through the integral. Hence, the mean square 
defiection contains the total number of states effecting on the electron level in the QD. The integral is limited above 
by kmax ~ In 3D-case it brings the factor R~^ ~ claimed in the beginning of the section. 

This result will change for other confined electric modes because of their different densities of states. This produces 
the different TV— scaling factor for the LS from these modes. 

The prefactor of the deflection, for any mode considered here, depends equally on the square root of the density of 
electrons, which is useful to be converted to rg, a characteristic length via the following deflnition: 27rrgn3D/3 = 1. 
Then, for 3D plasmon the deflection reads as: 

where the atomic length unit, as = /me^ ~ 0.53A, or the Bohr radius, gives the scale of the deflection (note that 
this definition docs not include any permittivity unlike an exciton Bohr radius in semiconductors). 

The depolarisation (the ratio of the level shift, 5E^ to the bare energy, E'^°'>) due to the 3D modes is as follows: 

Eio) 72^/67rV r, \R) ^ ' 
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The rude estimation of the prefactor shows that even for the small QD with = 100 the shift is 10~^ of the bare 
energy and will not be resolved because of a number of other different factors effecting the level position. 

To give a complete picture we note that the standard LS due to the zero-point oscillations of the free electromagnetic 
modes of the vacuum can be written as: 

^vac — I I 

TT \rs J 

where a ~ 1/137 is the fine structure constant, and the simple check shows that the logarithmic dependence of the 
last term in square brackets on N does not add any extra to the result and has to be dropped in our case. Though 
the slope of the LS in TV is much slower than in Eq.(||) the prefactor is tiny (~ 10"'') because of a^. 

IV. DEPOLARISATION: CONFINED MODES 

Let us consider the specific behavior of the LS materialized by the zero-point oscillations of the confined plasmon 
modes. The depolarisation in carbon shell cluster was shown jl] to be independent of the cluster size. The mean 
square deflection scales also as a zero power of the size 5'^ cx N. While it is interesting by itself, the carbon cluster 
matter will not be considered in the paper. However, there are confined modes in our QD problem those enhance the 
electrodynamical correction to the electron energy. 

With the decrease of the dimension of the field the plasmon density of states increases. Hence, the mean interaction 
of the electron with the plasmon field increases that will be evident from the scaling of (5^. 

Two possible candidates for the confined plasmon modes in the QD system, those have different densities of states, 
are the 2D plasmon and the OD spherical mode. The former mode can arise because of some interface possibly grown 
within the structure (see inset in the Figj^). It might be a conducting wetting layer, if it is thick enough to confine 
the electromagnetic field. The 2D plasmon naturally originates at the interface between the semiconductor structure 
and a metal . At the boundary of two dielectrics a surface plasmon is known to propagate . Its contribution will 
be discussed elsewhere as being smaller than 2D-plasmon one by a factor 10~^ at least owing to the fast space decay. 

The OD mode is the inherent property of the spherical inclusion of the foreign material in any matrix. The 
calculation of the frequency of this mode is slightly cumbersome (see Appendix for details). The surface QD mode 
has the quantum numbers L, M, the angular momentum and its projection on an axis, instead of 2D wave vector, k, 
for the standard 2D plasmon modes. The depolarisation is anomalous large in the OD case. It will be seen in this 
section from the scaling of S"^ and A. 

A. 2D plasmon 



The frequency of 2D plasmon is well known Q to depend on its 2D wave vector as: ujk = 2'Ke'^n2Dk / m. We 
will rewrite the 2D electron density, as before, in terms of the characteristic length: TTrln^D/l — 1 and perform the 
integration over the plasmon states. Then the mean square deflection can be expressed as: 




The scaling in N has a lower exponent that reflects the different density of the confined field (plasmon) states. 
Substituting 5"^ into the Hamiltonian given by Eq.(||), one gets the depolarisation as follows: 

The shift depends on the inverse size nearly linearly. However, the prefactor dominates at some moderate size of 
the QD and lessens the LS to 10""^ for N = 100. The depolarisation is still to be too small to expect experimental 
consequences. To be precise the result also depends on w, the distance between the 2D electrons and the QD. It is 
simply included in the consideration by multipling Eq.(|8|) by a factor ^/T:'Eivi{^J w / R) / {2^J w / R) ^ and the depolarisation 
declines 4 times at w/R^ 10. 



In 
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B. QD confined plasmon: Mode of cavity 



The b"^ considered above the less, the larger the QD size, that is not the case [Q for the giant deflection due 
to the completely localized modes. The localized modes are the surface plasmons of the spherical inclusion (with 
the dielectric function ei) in the matrix (with the different dielectric function 62). The frequency of the mode, wl, 
that we consider, is nearly the frequency of the bulk plasmon in the matrix, 10^2-, with the weak dependence on 
the mode angular momentum (see Appendix). The electric field of the zero-point oscillation is given by the formula 
■K{L + \/2)'hiOL/ . The summation over all states below some critical value Lc gives the mean square deflection: 
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(9) 



where it is natural to limit the summation above the excitation which wavelength is about the lattice constant d. We 
found PI that the 5'^ does not depend on the QD size: 



Sequently, the level shift depends on the size as R^^ (which comes from Eq.(|)): 

3^/6 V Ts \ d) \RJ 



(10) 



(11) 



Our estimation shows that the level correction, becoming of the order of 50%, plays the important role for the QD of 
100 atoms and smaller. We collected all studied contributions to the depolarisation LS and plot them in the log-log 
scale versus the QD size in Figure ^ 




FIG. 1. The level shift of the QD electron, calculated for 4 different depolarisation mechanisms. The giant shift, which 
results from the zero-point oscillations of the electric field of the QD surface mode, is shown in full line. The shifts by the 2D 
plasmon, the matrix bulk plasmon modes and the free field are shown as dash lines. The slopes and the prefactors of these 
depolarisation shifts are derived analytically and explained in the text. Inset: The scheme of the QD structure to model. 



The depolarisation because of the localized surface QD modes is large enough to propose an experiment supporting 
our model. It is easy to see that (5^ ~ w^'^, whence the LS depends on the mode frequency as well. Therefore, 
changing the optical properties of the matrix surrounding the QD, one shifts the levels. If the bare energy level, E'^°\ 
lies deep in the potential well, its position is nearly independent of the well depth which changes along with the matrix 
parameters. The deep bare level energy depends only on the well width R. Hence, keeping the same size of the 
QD and covering it with the different materials, one will derive solely the depolarisation LS, since it is distinguishable 
from the standard space quantization LS. 



V. SUMMARY 



The effect of the zero-point oscillations of the free and confined electromagnetic field on the level of the confined 
electron in the spherical QD is reviewed. The depolarisation due to an interaction with the zero-point oscillations of 
the field (produced by all other valence electrons) shifts up the bare one-electron state that seems to be a counterpart 
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for the vertex correction (electron-hole interaction, for example) which lowers the transition frequency down. It 
indicates that the studied effect should be taken into account for a many-body computation of a QD spectrum. 

To the best of our knowledge, the scaling dependence of the depolarisation level shift for the QDs is calculated in 
the first time. The size dependence of the LS is different for 4 cases considered in the paper. This scaling reflects 
that the different densities of states work in different mechanisms of the depolarisation due to the different 3D, 2D 
and OD-dimensional modes of the electric field are involved. Our model allows a theorist to skip a tedious quantum 
electrodynamical calculation but obtain the analytical selfconsistent estimation for the (many-body) level shift in a 
nanoscale system with the strong quantization. The result has not only a theoretical importance. 

Although, the depolarisation decreases with the QD size in general, the localized surface electromagnetic mode 
(which is specific to the QD as a void in the matrix material) results in the giant level shift and is to be possibly 
resolved experimentally for the QD made from some hundred atoms. Another method to detect the effect could 
be the measurement of a deep level position in the similar QDs buried by the substances with the distinct optical 
characteristics. Then the localized plasmon frequency changes along with the prefactor of the depolarisation shift, 
which could be observed by the optical spectroscopy of the QD system. 
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FIG. 2. The graphical solution for the QD plasmon modes with L > 1. The solid curve gives the function to solve (see 
secular equation in the text). The light dotted curves represent the dielectric functions of the dot and matrix substances. The 
full circles are the solutions for the sphere-like modes; the open circles are the cavity-like modes (see the text). Inset: The 
scheme of the band structure of a QD system. 



APPENDIX A: SURFACE QD PLASMON MODES 

The sought-for modes are given in complete spherical harmonics Pl{''')Yl.m{^), where Pl{i') is the Legendre 
polynomial and is the spherical harmonic The electrodynamic solution for the modes of the system consisting 
of the spherical particle with the dielectric function ei(w) and the surrounding matrix with the different dielectric 
function €2(0;) is one of the roots of the secular equation: 

44 = -^- (Al) 

ei(uj) L+l ^ ' 

The similar equation gives a mode of an empty void in the matrix in the limit ei = 1 (in the limit £2 = 1 it gives a 
mode of a sphere in the vacuum). The right hand side varies from —1/2 to —1. 

Let us suppose the simplest form for the dielectric function defined by one pole at some frequency, which is about 
the optical gap Eg, and one zero, which is given by the b ulk p lasma frequency of the material (both parameters differ 
for the materials 1 and 2). Then the solutions of the Eq. ([Al[ ) are easy to find graphically (see. Fig. 2). The function 
of w standing in the left hand of the expression is plotted in solid. It has two poles at Eg2 and ujpi and two zeroes at 
Egi and L0p2- All roots with L = 1, 2..Lc lie in between these two pairs as shown in the figure. The dependence of the 
mode frequency on the angular momentum is very weak. 

The lower-frequency modes, shown as black circles, correspond to the sphere-like plasmons. As lying above the 
QD gap, the lower mode effectively damps inside the QD. Therefore we will consider only the upper branch, shown as 
open circles, which is similar to the modes of the cavity. The frequency of the upper mode is close to the frequency 
of the bulk plasmon uip2 as it is seen from Figure (|^) . The mode frequency dependence on its angular momentum is 
negligible. 

The change of any of the 4 parameters, defining ei and £2, results in the QD mode frequency shift. However, for 
the sought-for cavity-like mode, the most important are the plasma frequencies of the materials. This provides the 
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mechanism for an experimental observation of the described effect in the matrix materials with the different cUp. It 
influences on the depolarisation LS, while the space quantization of the one-electron level, which is deep enough, 
depends solely on the QD width which has to be kept. 
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